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Abstract 

We calculate excitation functions for various dynamical quantities as well as exper- 
imental observables from SIS to RHIC energies within the HSD transport approach 
which is based on string, quark, diquark {q, q, qq, qq) and hadronic degrees of freedom 
without including any explicit phase transition to a quark-gluon plasma (QGP). It is 
argued that the failure of this more 'conventional' approach in comparison to exper- 
imental data should indicate the presence of a different phase which might be either 
attributed to space-time regions of vanishing scalar quark condensate {< qq > = 0) or 
to the presence of a QGP phase with strongly interacting partons. We study the K/ir 
ratio, the low mass dilepton enhancement in the invariant mass regime from 0.2 - 1.2 
GeV as well as charmonium suppression for central Au + Au collisions as a function of 
the bombarding energy and present predictions for these observables as well as hadron 
rapiditiy distributions at RHIC energies. Whereas all observables studied within HSD 
show smooth increasing/decreasing excitation functions, the experimental /tt^ ra- 
tio indicates a maximum at 11 A-GeV (or above) which is interpreted as a signature 
for a chirally restored phase in the course of the reaction. 
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1 Introduction 



The ultimate goal of relativistic nucleus-nucleus collisions is to reanalyze the early 'big-bang' 
under laboratory conditions and to find the 'smoking gun' for a phase transition from the 
initial quark-gluon plasma (QGP) to a phase characterized by an interacting hadron gas 

0, ^ . Any theoretical approach might describe such a phase transition starting from the 
partonic side with strongly interacting quarks and gluons ^ or from the hadronic side by 
involving hadronic degrees of freedom |]^, 0], i.e. hadrons with proper self-energies or spectral 
functions at high baryon density or temperature. It remains to be seen which approach will 
prove to be more successful, economic and transparent. 

Nucleus-nucleus collisions with initial energies per nucleon of ^/s = 200 GeV or ^21.5 A-TeV 
will be available soon at the Relativistic-Heavy-Ion-Collider (RHIC) in Brookhaven. In cen- 
tral collisions of Au -|- Au here energy densities above 5 GeV/fm'^ are expected such that 
the critical energy density for a QGP phase should be overcome in considerable space-time 
volumes where the relevant degrees of freedom are partons (quarks and gluons). Parton cas- 
cade calculations have been used so far 0, ^ to estimate the energy densities and particle 
production yields in violent reactions at y/s = 200 GeV, an order of magnitude higher than at 
SPS energies (a/s ~ 20 GeV). Intuitively one expects that the initial nonequilibrium phase 
of a nucleus-nucleus collision at RHIC energies should be described by parton degrees of 
freedom whereas hadrons are only formed (by 'condensation') at a later stage of the reaction 
which might be a couple of fm/c from the initial contact of the heavy ions. Thus parton cas- 
cade calculations - including transitions rates from perturbative QCD - should be adequate 
for all initial reactions involving a large 4-momentum transfer between the consituents since 
QCD is well tested in its short distance properties. The question, however, remains to which 
extent the parton calculations can be extrapolated to low where hadronic scales become 
important. As a rough estimate one can employ here the average mass of vector mesons, the 
nucleon and its first excited state, which gives Q^rit ~ ^ GeV^. On the other hand, using the 
uncertainty relation this implies time scales of At < 0.2 fm/c = Qcrit relative separations 
of partons Ar < 0.2 fm, which are small compared to the hadronic size or average life time 
of the p, A, etc. in free space or the formation time of hadrons tp ^ 0.7 — 0.8 fm/c as used 
in the HSD transport approach |TD|, [TT], |T^. 

Turning the argument around, a nonequilibrium hadronic approach involving a time scale 
of 0.7-0.8 fm/c cannot tell anything about shorter times because the uncertainty relation does 
not allow to distinguish states which are separated in mass by less than ^ 300 MeV = tp^, 
which is the N — A mass difference. Thus one faces the problem that neither the parton 
description nor a nonequilibrium hadronic model should be valid for times 0.2 fm/c < t < 
0.7-0.8 fm/c in individual hadronic reactions, which corresponds to the nonperturbative 
formation time of the hadronic wavefunction. This regime of the 'soft' QCD physics is 
presently not well understood and appropriate dynamical models are urgently needed. In 
the HSD approach these intermediate times are described by color neutral strings, where the 
leading quarks and 'diquarks' in a baryonic string (or quarks and antiquarks in a mesonic 
string etc.) are allowed to rescatter again with hadronic cross sections divided by the number 
of constituent quarks and antiquarks in the hadrons, respectively [p!^ . 

Furthermore, the question of chiral symmetry restoration at high baryon density and/or 
high temperature is of fundamental interest, too 0. Whereas lattice QCD calculations at 
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zero baryon chemical potential indicate that a restoration of chiral symmetry goes along with 
the deconfinement phase transition at some critical temperature T^, the situation is less clear 
for finite baryon density where QCD sum rule studies indicate a linear decrease of the scalar 
quark condensate < qq > - which is nonvanishing in the vacuum due to a spontaneous 
breaking of chiral symmetry - with baryon density ps towards a chiral symmetric phase 
characterized by < gg > = 0. This decrease of the scalar condensate is expected to lead to 
a change of the hadron properties with density and temperature, i.e. in a chirally restored 
phase the hadrons might become approximately massless |jl3| or at least vector and axial 
vector currents should become equal 0; the latter implies that the p and ai spectral 
functions should become identical. Since the scalar condensate < gg > is not a direct 
observable, its manifestations should be found in different hadronic abundancies and spectra. 

Nowadays, our knowledge about the hadron properties at high temperature or baryon 
density is based on heavy-ion experiments from BEVALAC/SIS to SPS energies where hot 
and dense nuclear systems are produced on a timescale of a few fm/c. As mentioned above, 
the information from ultrarelativistic nucleus-nucleus collisions at RHIC, i.e. initial y/s = 200 
GeV per nucleon, will be available soon . However, any conclusions about the properties of 
hadrons in the nuclear environment are based on the comparison of experimental data with 
nonequilibrium kinetic transport theory |]TU|, |TB|, |TB|, |T^. As a genuine feature of transport 
theories there are two essential ingredients: i.e. the baryon (and meson) scalar and vector 
self-energies as well as in-medium elastic and inelastic cross sections for all hadrons involved. 
Whereas in the low-energy regime these 'transport coefficients' can be calculated in the Dirac- 



Brueckner approach starting from the bare nucleon-nucleon interaction [g^, 21, 22| this is no 



longer possible at high baryon density (ps > 2-3po) and high temperature, since the number 
of independent hadronic degrees of freedom increases drastically and the interacting hadronic 
system should enter a phase with < gg >^ [jl3|, ^ ^ ^ as mentioned before. As a 
consequence the hadron self-energies or spectral functions in the nuclear medium will change 
substantially especially close to the chiral phase transition and transport theoretical studies 
should include the generic properties of QCD in line with nonperturbative computations on 
the lattice ||, |^ ||, ||]. 



In this work we concentrate on excitation functions of hadronic observables from SIS 
to RHIC energies with the aim to find out optimal experimental conditions to search for 
'traditional' phenomena such as strangeness enhancement (Section 3), low mass dilepton 
enhancement (Section 4) or charmonium suppression in nucleus-nucleus collisions (Section 
5). Our studies are performed within the HSD transport approach that has been described 



in Refs. 12 1 and been tested for p -\- p, p -\- A and A + A collisions from the SIS to SPS 



energy regime [0. Actual predictions for hadron rapidity spectra and J/\E' suppression will 
be presented in Section 6 while Section 7 concludes the study with a summary. 



2 Theoretical considerations 

In this Section we briefiy recall the ingredients of the covariant transport theory, that has 
been denoted as Hadron-String-Dynamics (HSD) [|1^, which formally can be written as 
a coupled set of transport equations for the phase-space distributions fh{x,p) of hadron h 
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m, 11, m, i.e. 



= /rf2d3rf4...[G'tG']i2^34...5r(n + n2-n3-n4...) 

h2h2.hi... 

X {fh3{x,P3)fh^{x,P4:)fh{x,p)fh2{x,P2) 

- fh{x,p)fh2{x,P2)fh;{x,P3)fh^{x,P4)} ... . 



In Eq. (|l|) U^{x,p) and U\^{x^p) denote the real part of the scalar and vector hadron self- 
energies, respectively, while [G''''G]i2-»34...'^r(n + 112 — Ha — 114 . . .) is the 'transition rate' for 
the process 1 + 2— *-3-|-4-|-... which is taken to be on-shell in the semiclassical limit adopted. 
The hadron quasi-particle properties in ([I|) are defined via the mass-shell constraint , 

5{Ii,W-Mf) , (2) 

with effective masses and momenta (for a hadron of bare mass and momentum p^) given 
by 

Ml{x,p) = Mh + U^{x,p) 

U^{x,p) = p^-U}:{x,p) , (3) 

while the phase-space factors 

fh{x,p) = l±Mx,p) (4) 

are responsible for fermion Pauli-blocking or Bose enhancement, respectively, depending 
on the type of hadron in the final/initial channel. The dots in Eq. (0) stand for further 
contributions to the collision term with more than two hadrons in the final/initial channels. 
The transport approach (|1|) is fully specified by U^{x,p) and Uf^{x,p) (/x = 0, 1, 2, 3), which 
determine the mean-field propagation of the hadrons, and by the transition rates G^G S^{. . .) 
in the collision term, that describe the scattering and hadron production/absorption rates. 
In Ref. |jll[ the scalar and vector mean-fields and Uj^ have been determined in the 



mean-field limit from an effective hadronic Lagrangian density Ch that has been fitted to 
the equation of state of nucleonic matter as resulting from the Nambu-Jona-Lasinio (NJL) 
model. Without going through the arguments again we show in Fig. |T| the nucleon scalar 
(Us) and negative vector potential (— f/o) as a function of the nuclear density p and relative 
momentum p of the nucleon with respect to the nuclear matter rest frame. Whereas the 
vector potential increases practically linearly with density (at low momentum p) the scalar 
potential saturates with density such that the nucleon effective mass M* = Mq + Us almost 
drops to zero for p > 0.6 fm~^. Both potentials decrease rather fast in magnitude with 
momentum p and practically vanish above a few GeV/c. 
In Fig. the real part of the potential 



UsEP = Uoipo,p) + ^p' + {MN + Us)'-^p' + Mj, (5) 

is shown again as a function of p and p. Whereas we see a net attraction for momenta |p| < 
0.5 GeV/c up to densities of ~ 0.3 fm~^, the net potential becomes repulsive for higher 
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momenta, reaches a maximum repulsion at |p| ~ 1 GeV/c and then drops again with |p|. 
We mention that at density po the potential Usep compares well with the potential from the 
data analysis of Hama et al. as well as Dirac-Brueckner computations from up to 



a kinetic energy E^in of 1 GeV |]rT|. The formula (^) reduces to the familiar expression for 



the Schroedinger equivalent potential (Eq. (3.16) of Ref. []IU[) in the low density limit. 

In our transport calculations we include nucleons, A's, N*(1440), N*(1535), A, E and 
E* hyperons, S's and f2's as well as their antiparticles. In a first approximation we assume 
here that all baryons (made out of light {u, d) quarks) have the same scalar and vector self- 
energies as the nucleons while the hyperons pick up a factor 2/3 according to the light quark 
content and S's a factor of 1/3, respectively. 

In the HSD approach the high energy inelastic hadron-hadron collisions are described by 



the FRITIOF model |33], where two incoming hadrons emerge the reaction as two excited 



color singlet states, i.e. 'strings'. According to the Lund model [^] a string is characterized 
by the leading constituent quarks of the incoming hadron and a tube of color flux is supposed 
to be formed connecting the rapidly receding string-ends. In the HSD approach baryonic 
{qq — q) and mesonic (g — q) strings are considered with different flavors (g = m, d, s). In the 
uniform color field of the strings virtual qq or qqqq pairs are produced causing the tube to 
fission and thus to create mesons or baryon-antibaryon pairs. The production probability P 
of massive ss or qqqq pairs is suppressed in comparison to light flavor production (mm, dd) 
according to a Schwinger-like formula |3^ 



— — = 7, = exp -TT ^, 6 

F[uu) \ Ik ) 

with K «i 1 GeV/fm denoting the string tension. Thus in the Lund string picture the 
production of strangeness and baryon-antibaryon pairs is controlled by the constituent quark 
and diquark masses. Inserting the constituent quark masses ra^ = 0.3 GeV and = 0.5 
GeV a value of 7s ~ 0.3 is obtained. While the strangeness production in proton-proton 
collisions at SPS energies is reasonably well reproduced with this value, the strangeness yield 
for p + Be collisions at AGS energies is underestimated by roughly 30% (cf. ||12[)- For that 
reason the relative factors used in the HSD model are |T^ 

_ J 1 : 1 : 0.3 : 0.07 , at SPS to RHIC energies , . 

~ \ 1 : 1 : 0.4 : 0.07 , at AGS energies, ^ ' 

with a linear transition of the strangeness suppression factor 7^ = s : ti as a function of y/s 
in between. 

Additionally a fragmentation function f{x, rrit) has to be specified, which is the prob- 
ability distribution for hadrons with transverse mass rrit to acquire the energy-momentum 
fraction x from the fragmenting string. 



f(x, rrit) ~ —(1 — exp (—bm'f/. 
X ^ 

with a = 0.23 and b = 0.34 GeV^^ p|. 



We recall that the LUND model includes partonic diffractive scattering and mini- 



jet production as well [^. The latter phenomena are not important at SPS energies and 
below, however, become appreciable at RHIC energies. In this respect the HSD approach 
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dynamically also includes the hard partonic processes as far as quarks and antiquarks are 
involved. However, it does not employ hard gluon-gluon processes beyond the level of 'string 
phenomenology'. This has to be kept in mind with respect to the predictive power of the 
model at RHIC energies and beyond. 

The medium modifications due to the hadron self-energies, furthermore, require to in- 
troduce some conserving approximations in the collision terms in line with the modified 
quasi-particle properties. Since these in-medium modifications - related to 'low momentum 
physics' - are not of primary interest in this study we discard an explicit discussion here and 



refer the reader to Ref. HTl 



2.1 The scalar condensate 

The scalar quark condensate < gg > is viewed as an order parameter for the restoration of 
chiral symmetry at high baryon density and temperature. A model independent relation for 
the scalar quark condensate at finite (but small) baryon density and temperature has been 
given by Drukarev and Levin i.e 



<qq> , ps 



<qq>v fn^nl 



(9) 



where < qq >v denotes the vacuum condensate, ^ 45 MeV is the pion-nucleon S-term, 
/tt and m^r the pion decay constant and pion mass, respectively, while E{p)/A is the binding 
energy per nucleon. Since for low densities the scalar density ps in @ may be replaced by 
the baryon density ps, the change in the scalar quark condensate starts linearly with ps and 
is reduced by a factor 1/3 at saturation density po- A simple linear extrapolation then would 
indicate that at ~ 3po a restoration of chiral symmetry might be achieved in heavy-ion 
collisions. 

A reasonable estimate for the scalar condensate in dynamical calculations has been sug- 
gested by Friman et al. [P7[, 



<qq>v fn^i V fn^i 

where ah denotes the cr-commutator of the relevant mesons h. For pions and mesons made 
out of light quarks and antiquarks we use ah = whereas for mesons with a strange 

(antistrange) quark we adopt ah = ttIt^/A according to the light quark content. Within the 
same spirit the E-term for hyperons is taken as 2/3 ~ 30 MeV while for S's we use 1/3 
^ 15 MeV. 

The scalar density of mesons (of type h) is given by 



, (2g + l)(2t + l) f^, rrih , , 



+ m 



with fh denoting the meson phase-space distribution of species h. In (|n|) s, t denote the spin 
and isospin quantum numbers, respectively. We note that the scalar density of baryons is 
calculated in line with ( [TT| ) by replacing the mass and momentum by the effective quantities 
in m. 
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The actual numerical result for the space-time dependence of the scalar condensate ([T0|) 
is shown in Figs. ^ and | for central Au + Au collisions at 6 A-GeV and 20 A-GeV, re- 
spectively. Here the condensate is divided by the vacuum condensate < qq >v such that 
the nonperturbative vacuum is characterized by a value of 1. Furthermore, the 2;-direction 
has been stretched by the Lorentz-factors to compensate for Lorentz contraction, while 
negative numerical values for the condensate have been suppressed. According to (H), ([T0|) 
the scalar condensate is reduced inside the approaching nuclei by about 35%; this reduction 
becomes more pronounced when the nuclei achieve full overlap. As seen from Fig. ^ al- 
ready at 6 A-GeV there is a substantial space-time region of vanishing scalar condensate for 
6 fm/c < t < 13 fm/c, where the conventional hadronic picture is not expected to hold any- 
more. The space-time volume of vanishing quark condensate slightly increases for 20 A-GeV 
(Fig. ^), however, the increase is only moderate and resembles very much the situation for 
the space-time integral of high density baryon matter (cf. Fig. 1 in Ref. |]38|). 

We mention that at higher bombarding energies the 4- volume {x = {t, r)) 

V{a) = [ dhdt e{a - (12) 

that counts the fraction of the scalar condensate below the value of a (0 < a < 0.3) is 
essentially determined by the pion density whereas below about 20 A-GeV the dominant 
contribution stems from the baryons (cf. also Ref. ||3^). Since the pion density can be 
considered as a measure of the vacuum 'temperature', a chiral order transition below about 
20 A-GeV is dominated by the baryon density while especially at SPS or even RHIC energies 
a chiral order transition is due to 'temperature'. 

The question now arises, if there are proper experimental observables that allow to trace 
down such type of phase transition (or cross over). When gating on central collisions of 
Au -|- Au (or Pb + Pb) such phenomena should show up in their excitation functions. We 
note that in a pure hadronic transport approach we expect a smooth behaviour of practically 
all observables with bombarding energy due to an increase of thermal excitation energy. This 
is not so obvious for the HSD approach where a gradual transition from hadronic excitations 



to strings and quark (or diquark) degrees of freedom is involved. In fact, in Ref. it 
has been claimed that the excitation function of transverse and elliptic proton flow together 
with the transverse pt spectra of protons suggest a transition from hadron to string matter 
at about 5 A-GeV. Here, however, we concentrate on meson abundancies and spectra and 



refer the reader to Ref. for the collective dynamical aspects and to Ref. |^ for the 
thermal properties of the theory that involves a limiting ('Hagedorn') temperature of Ts ~ 
150 MeV due to the continuum string excitations. 



3 Meson production 

To present a general overview on the meson abundancies in nucleus-nucleus collisions we 
show in Fig. ^ the meson multiplicities for central collisions of Au + Au from SIS to RHIC 
energies. All multiplicities for tt'^ ,r], , K~ ,(f) as well as J/\l/ mesons show a monotonic 
increase with bombarding energy which is only very steep at 'subthreshold' energies, i.e. at 
bombarding energies per nucleon below the threshold in free space for NN collisions. At 
higher bombarding energies the meson abundancies group according to their quark content. 
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i.e. the multiplicities are reduced (relative to tt"*") by about a factor of 4-5 for a strange 
quark, a factor of ~ 50 for ss = and a factor of ~ 2 ■ 10^ for cc = J/\l/. We mention that in 
these calculations the meson rescattering and reabsorption processes have been taken into 
account; this reduces the J/\l/ cross section at RHIC energies by about a factor of 10 (cf. 
Section 5). At 'subthreshold' energies the multiplicity turns out to be almost the same as 
the antikaon multiplicity, but then rises less steeply with bombarding energy. Apart from 
the excitation function - that still has to be controlled experimentally - we thus find no 
pronounced change in the shape of the meson abundancies up to RHIC energies where data 
are expected to come up soon. 



We recall that our investigations on strangeness production up to 2 A-GeV |41, H2l have 



given some evidence for attractive antikaon potentials in the medium while for kaons only 



a very moderate repulsive potential was suggested |H2|; rj mesons apparently do not show 



sensible in-medium effects according to the studies in Ref. Q in comparison to the 
available experimental spectra. At AGS and SPS energies, on the other hand, the potential 
effects on kaon and antikaon abundancies have been found to be only very low [jlO[ such that 
meson potentials have been discarded in Fig. |^ for the overview. 

Since the meson abundancies show no sudden change in the excitation function, we turn 
to particle ratios. Here strangeness enhancement has been suggested for more than 2 decades 
to possibly qualify as a sensible probe for a QGP phase (cf. Ref. [0 for a recent overview). 
Here we examine the K^/n^ ratio at midrapidity in central Au + Au (or Pb + Pb) collisions 
where experimental data are now available from SIS to SPS energies [^, |^, We recall 
that detailed comparisons of pion and kaon rapidity distributions and transverse momentum 



spectra to the available data have been presented in Refs. [jTO|, such that we directly can 
continue with the corresponding excitation function. 

In order to discuss the strangeness production over the complete energy range from SIS to 
RHIC energies we show in Fig. ^ the calculated K'^/tt^ ratio at midrapidity (-0.5 < i/cm < 
0.5) for central Au + Au collisions in comparison to the experimental data. This ratio 
experimentally is substantially lower at SPS energies (~ 13.5%) compared to AGS energies 
(~ 19%). At SPS energies this ratio is only enhanced by a factor 1.75 for central Pb + Pb 
collisions relative to p + p reactions and has to be compared to the factor ^ 3 at AGS. 
Such a decrease of the scaled kaon yield from AGS to SPS energies is not described by 
the HSD transport model (without kaon self-energies) which shows a monotonic increase 
with bombarding energy similar to pp reactions (open circles). Furthermore, the higher 
temperatures and particle densities at SPS energies tend to enhance the /n^ yield closer 
to its thermal equilibrium value of ~ 20 — 25% |^ at chemical freezout and temperatures 
of T ^ 150 MeV. 

Our findings have to be compared to results obtained by other microscopic approaches. 
Here only the RQMD model very recently provides a study partly comparable to that 
of Ref. [0], however, excluding p -|- A reactions to control the amount of production by 
rescattering. Whereas kaon production in p + p reactions is comparable to our results in [|l^ 
(input of the transport model), the RQMD model yields a higher jn^ ratio in Au + Au, 
Pb + Pb collisions at all energies due to a higher kaon yield from rescattering. The latter 
can be attributed to 'high mass strange resonances' that have been incorporated to describe 



the low energy kaon production via resonance production and decay |Q (s-channels). Since 



these 'high mass resonances' can be repopulated in resonance-resonance scattering, a rather 
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high ^/s is concentrated in a single degree of freedom for a short time. Such 'hot spots' 
then lead to a higher /tt^ ratio in A + A reactions especially at AGS energies. While 
the /tt^ ratio can be reasonably described at AGS energies in Au + Au reactions, this 
ratio is overestimated significantly at SPS energies (cf. Fig. 4 of Ref. |^). Thus also the 
RQMD model does not describe the relative decrease of the j-n^ ratio from 11 A-GeV 
to 160 A-GeV. The results of the RQMD calculations at SIS energies are not known to the 
authors. 



4 Dilepton production 



Electromagnetic decays to virtual photons (e^e or pairs) have been suggested long 

ago to serve as a possible signature for a phase transition to the QGP [0, |4|, ^ 



or to be an ideal probe for vector meson spectroscopy in the nuclear medium 
out in Refs. 
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As pointed 

Q the isovector current-current correlation function is proportional to the 
imaginary part of the p-meson propagator and also to the dilepton invariant mass spectra. 
Dileptons are particularly well suited for an investigation of the violent phases of a high- 
energy heavy-ion collision because they can leave the reaction volume essentially undistorted 
by final-state interactions. Indeed, the dilepton studies in heavy-ion collisions by the DLS 
Collaboration at the BEVALAC ^ and by the CERES |5|,|5|, HELIOS m P, NA38 [||] 
and NA50 Collaborations |Q at SPS energies have found a vivid interest in the nuclear 
physics community. 

We recall that the question of chiral symmetry restoration does not necessarily imply that 
vector meson masses have to drop with baryon density or temperature 0] . Actually, chiral 
symmetry restoration (ChSR) only demands that the isovector current-current correlation 
function and the axial vector current- current correlation function (dominated by the chiral 
partner of the p, the Oi-meson) should become identical at high pb or temperature T, respec- 
tively, because there should be no more differences between left- and right-handed particles 
or equivalently vector and axial vector currents |l^, Thus also a strong broadening of 
the p- as well as the Oi-spectral function and their mixing in the medium ||6^, |6^, Q can 
be considered as a signature for ChSR which, however, is not easy to detect experimentally. 



In Ref. 



we have demonstrated that the present experimental data on the low mass 
dilepton enhancement at SPS energies can be described equally well within the 'dropping 
p-mass' scenario as well as within the 'melting' p picture, which implies a large spreading 
in mass of the p spectral function due to its couplings to baryons and/or mesons. The 
situation at SIS/BEVALAC energies is more 'puzzling' since here the low mass dilepton 
enhancement is neither described in the dynamical spectral function approach nor in 
the 'dropping mass' scheme though the dilepton data from the DLS collaboration 
are reproduced within the known sources rather well from 1-5 GeV bombarding energy 
|]7ll , |72| . In short, the dynamical origin of the low mass dilepton enhancement is not yet 
understood. 

Here we propose to investigate the excitation function of low mass dilepton enhancement 
in central Au + Au collisions. As discussed in Refs. |63] this excess of dileptons is most 



probably due to the isovector (p) channel, i.e. the imaginary part of the isovector current- 
current correlation function which, however, mixes with the axial vector current-current 
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correlation function at finite temperature and baryon density |6^ 



Since the p-meson spectral function is of primary interest, it is important to have some 
information about the actual baryon densities that the p-meson experiences during its prop- 
agation and decay in central Au + Au collisions. This information is displayed in Fig. ^ 
- as resulting from the HSD transport calculation - for central reactions at 2, 5, 10, 20, 
40 and 160 A-GeV. Here the meson-baryon production channels and baryon-baryon produc- 
tion channels are summed up by the solid lines and denoted as [ttB —>■ p,BB p). The 
meson-meson production channels such as tttt — *■ p, Oi — >■ 7rp etc. are summed up in the 
dashed histograms indicated as tttt — p according to the dominating channel. We find that 
especially the initial BB production channels produce p-mesons at very high densities close 
to 2 po7cm where 7cm is the Lorentz factor in the cms. However, these production channels 
are much less abundant than the meson-meson channels which extend over a larger time 
span (for the heavy system Au + Au) and essentially occur at much lower baryon density, 
respectively. This effect becomes even more pronounced with increasing bombarding energy, 
i.e. 40 - 160 A-GeV, where most of the p-mesons are produced at baryon densities below 
2 pq. Since in Fig. |^ the relative abundancy of p-mesons is displayed as a function of the 
baryon density at the production (formation) point, the time averaged value of the density 
is even lower due to a fast expansion of the hadronic fireball. Thus to probe on average high 
baryon densities by p-mesons one should step down in energy to 2 - 5 A-GeV in order to 
optimize effects due to the coupling to baryons. 

A general overview of dilepton mass spectra in central (6=2 fm) Au + Au collisions from 
2 A-GeV to 21.5 A-TeV is presented in Fig. where the upper part corresponds to the 
case of vacuum spectral functions for all mesons (cf. Fig. 8.20 of Ref. [|1^), while the lower 
part is calculated by employing the p spectral function from Rapp et al. Whereas 



for the free meson spectral function one observes essentially an increase of the dilepton 
production channels with bombarding energy without any substantial change in the spectral 
shape (except for an increasing peak from the 0; cf. Fig. 5), the in-medium calculations 
yield almost exponential mass spectra above M > 0.4 GeV with small peaks from vacuum 
LJ and decays at M f« 0.78 GeV and 1.02 GeV. 

We mention that for the vacuum spectral functions (upper part of Fig. 8) the shape of the 
dilepton spectra (for Mg+g- > 0.15 GeV) is due to the superposition of rj, rj', uj and ai Dalitz 
decays and direct vector meson decays {p,u!,(f)), where all mesons may also be produced 
in meson-baryon and meson-meson collisions. The increase in dilepton yield (for Mg+e- > 
0.15 GeV) with bombarding energy thus is due to an enhanced production of rj, rj', p, uj, 0, oi 
mesons. Their relative abundance from SIS to RHIC energies does not scale directly with 
the charged particle multiplicity which is dominated by protons, vr^ and (cf. Fig. 5). 
However, above about 50 - 100 A-GeV the meson ratios do not change very much (cf. Figs. 
5,6) while the charged particle multiplicity becomes dominated by vr^ and K^. Thus from 
SPS to RHIC energies the low mass dilepton yield should approximately be proportional to 
the charged particle multiplicity, too. 

The relative change in the dilepton spectra is quantitatively displayed in Fig. |^ - again 
for central {b= 2 fm) collisions of Au -|- Au - for the case of free meson spectral functions 
(solid lines) and the spectral function from Rapp et al. (dashed lines). In line 

with the discussion above the most prominent spectral changes are expected at rather low 
bombarding energies from 2-10 A- GeV, where the enhancement in the invariant mass range 
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from 0.3-0.6 GeV is about a factor of 4. Note that the uj and vacuum decays show clear 
peaks on top of the 'exponential' continuum, which will have to be identified experimentally. 
We mention that in all our calculations we have incorporated an experimental ('optimistic') 
mass resolution of AM = 10 MeV for the dilepton invariant mass. 



5 Charmonium production and suppression 



Matsui and Satz |]73[ have proposed that a suppression of the J/\E' yield in ultra- relativistic 
heavy-ion collisions is a plausible signature for the formation of the quark-gluon plasma 
because the J/\l/ should dissolve in the QGP due to color screening. This suggestion has 
stimulated a number of heavy-ion experiments at CERN SPS to measure the J/\l/ via its 
dimuon decay. Indeed, these experiments have shown a significant reduction of the J/\l/ yield 
when going from proton- nucleus to nucleus- nucleus collisions |63]. Especially for Pb + Pb 



at 160 A-GeV an even more dramatic reduction of J/\l/ has been reported by the NA50 
Collaboration 



To interpret the experimental results, various models based on J/\E' absorption by hadrons 
have been also proposed (cf. Refs. |IU| , [T7| , [T^ for recent reviews) that do not involve the 
assumption of a QGP phase transition. The role of comover dissociation is presently again 
heavily debated [|, |78|] especially since theoretical calculations for J/\l/-meson dissociation 
cross sections differ by up to a factor of 50 0, |80|, ^ The problem is even more 
complicated since the J/\l/ meson is not created instantly as a hadronic state and there is 
also a substantial (~ 35 %) feeding from the Xc — 7 decays. Moreover, it is expected that the 
cc pair is first produced in a color-octet state together with a gluon ('pre- resonance state') 
and that this more extended configuration has a larger interaction cross section with baryons 
and mesons before the J/\l/ singlet state, or Xc finally emerges after some formation time 
Tc. Additionally, the dissociation on mesons of formed J/\l/'s will differ from Xc due to their 
different thresholds with respect to the DD channel as well as for the Since experimental 
information on the various charmonium-meson cross sections - especially at low relative 
momenta - will be hard to obtain, the excitation function of charmonium suppression in 
central nucleus-nucleus collisions might be exploited to obtain additional information on the 
absorption scenarios. One expects that quite below the bombarding energy necessary for the 
formation of a QGP phase the charmonium absorption should be entirely due to dissociation 
with hadrons; any additional suppression due to color screening then will show up in a more 



rapid suppression with the incident energy or with the energy density achieved |]83 |. 

In this Section we calculate the excitation function for J/\l/ suppression within two ab- 
sorption scenarios, i.e. the 'early' and 'late' comover dissociation models, which have been 
explored in detail by our group before at SPS energies ^ The 'early' comover 
absorption scenario is based on the idea, that a cc pair is created in the initial 'hard' phase 
of the nucleus-nucleus collision, where the string density is very high, and the cc pair is 
dissolved in the color electric field of neighboring strings. Since in the HSD approach the 
information on the string density as well as the string space-time extension is available, the 
absorption model has only a single parameter, i.e. the average transverse dimension of an 
extending string. In Ref. |86| a string radius of 0.2 fm was found to describe simultaneously 
the data for p + A and S + U at 200 A- GeV from NA38 P[ and for Pb + Pb at 160 A- GeV 
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from NA50 ||63[ when adopting the conventional charmonium-nucleon dissociation cross sec- 
tion of 6 mb. It has been also speculated that the overlap of strings due to percolation might 
describe the phase transition to a QGP [p^ ]. 

In the 'late' comover scenario the additional charmonium suppression is due to charmo- 
nium-meson scattering to DD with an average charmonium-meson cross section of ~ 3 mb 
| 85| . Cross sections of this order have been calculated by Haglin in Ref. |j79| within a 
meson-exchange model and thus might appear not unrealistic. In our present calculation we 
refer to the model II of Ref. including a 'pre-resonance' charmonium life time of 0.3- 
0.5 fm/c which is supported (within the errorbars) by a more recent analysis of charmonium 
suppression as a function of the Feynman variable Xi? in p + A reactions |^8[ from He et 
al. and Kharzeev et al. |Q. We recall that within the model II of Ref. [^] the J/'^ 
suppression data for p + A and S -|- U at 200 A-GeV from NA38 [^| and for Pb -|- Pb at 
160 A-GeV from NA50 ^ 



63 



have been described very well when adopting a 'pre-resonance '- 
nucleon cross section of 6 mb and a J/\E'-nucleon cross section of 3-4 mb in line with the data 
on J/\E' photoproduction |91|. Independent dynamical studies on charmonium suppression 
within the UrQMD model later on have lead to very similar conclusions. 

Since the comparison of our calculations at SPS energies has been performed to data 

TO a comparison of the 'early' comover model 
with the more recent data from NA50 

IB 



taken in 1995 and before we show in Fig. 
(dashed line) and the 'late' comover model II ||85[ 

|]76|, pi] for Pb + Pb at 160 A-GeV using the explicit numbers from Refs 



85 



The 

'early' comover absorption model here gives a little too low suppression at high Et whereas 
the 'late' comover absorption model is still in line with the more recent data from 1996 
and 1998 with minimum bias (open triangles and open circles). The data from 1996 (full 
squares), that show a (much debated) two-step behaviour, do not agree with the explicit Et 
dependence from our calculations; however, the 1996 minimum bias data (open triangles) 
well match for Et > 40 GeV whereas the J/^ suppression is slightly overestimated at lower 
Et- We do not comment on the highest Et data points from 1998 with minimum bias since 
our earlier analysis did not extend to these specific events. 

The question now arises, if the excitation function for J/\E' suppression in central Au + Au 
collisions might show some unusual behaviour within the two scenarios discussed above or 
how they might be disentangled. In order to achieve the same suppression factor in central 
collisions within the 'early' comover model we have increased the string absorption radius 
from rs = 0.2 fm to 0.22 fm to get the same value for the J/^f survival factor Sj/^ at 
160 A- GeV. We note that the total J/"^ multiplicity shown in Fig. ^ has been calculated 
within the 'late' comover model. It drops by almost 3 orders of magnitude when decreasing 
the bombarding energy from 160 A-GeV to 20 A-GeV. At the lowest energy considered 
here the experimental J/"^ signal will be very hard to measure; it is hopeless within the 
present experimental setups. Nevertheless, it is worth exploring theoretically if some unusual 
excitation function might be found. Our results are displayed in Fig. [Tl| (l.h.s.) for the 'early' 
absorption model (open circles) and for the 'late' comover model (full circles); both models 
practically do not differ in their excitation functions and show a very smooth decrease of 
Sj/ij, from 0.4 to 0.3 with increasing bombarding energy. The net absorption by baryons is 
dominant in both scenarios, however, differs in magnitude due to the simple fact that in the 



^This comparison is necessary since the 1995 data have been rescaled in |94| and our calculations are 
reported inconsistently in the more recent presentations of this topic |94, M . 



12 



'early' (string) absorption model there is less suppression by baryons since the absorption 
by strings competes at early times. In the 'late' comover model there is more absorption 
by baryons because the mesons are formed at later stages and not competitive in the early 
phase; their relative contribution is lower as for strings accordingly. However, these individual 
contributions cannot be distinguished experimentally and thus are 'irrelevant'. 

The r.h.s. of Fig. 11 shows the survival probability Sj/^i in the 'late' comover model for 
central Au + Au collisions from 0.160 A-TeV to 21.5 A-TeV, respectively, where we have 
gated on J/^^'s in the rapidity interval -1 < Ucm < 1- The solid line stands for the total J /'^ 
survival probability while the dashed line displays the relative absorption on baryons and 
the dotted line the relative dissociation on mesons. Whereas the dissociation on baryons 
is practically constant with bombarding energy, the absorption on mesons increases in line 
with the higher meson densities achieved with increasing E/A. We note that the 'early' 
comover model leads to a similar total absorption within the numerical accuracy. 

We have to mention that neither the 'early' nor the 'late' comover model might be 
realized in nature exclusively and both absorption processes should occur within the same 
reaction with probably different weights. Since we do not find a substantial difference for 
both scenarios also a linear combination of both absorption models, i.e. decreasing as 
well as the charmonium-meson cross section accordingly, will lead to a similar excitation 
function. This also holds for the relative suppression on the transverse energy Ex (cf. Fig. 

Inspite of the rather disappointing perspectives to disentangle the 'late' and 'early' co- 
mover models experimentally at the full range of SPS energies, the excitation function of 
J/^P still might show some discontinuity m E/A experimentally, which could rule out the 
two models studied here and indicate a transition to a QGP phase. This subject is taken up 
in the next Section again with respect to the dependence of Sj/^, on the transverse energy 
Et produced in Au + Au colhsions at RHIC energies. 



6 Predictions for RHIC energies 

As noted in the introduction one expects that the initial nonequilibrium phase of a nucleus- 
nucleus collision at RHIC energies should be described by parton degrees of freedom, whereas 
hadrons are only formed (by 'condensation') at a later stage of the reaction and interact 
until freeze out. Thus parton cascade calculations should be adequate for all initial reactions 
involving a large 4-momentum transfer between the constituents, while hadron cascades 
should be appropriate in the final hadronic expansion phase. We suggest that the dynamics 
in between the partonic and hadronic phase might be described by quarks (diquarks) and 
strings as e.g. implemented in the HSD approach. 

The practical question is, however, if nonequilibrium partonic and hadron/string mod- 
els can be distinguished at all, i.e. do they lead to different predictions for experimental 
observables? In fact, first applications of the parton cascade model developed by Geiger 
P, to nucleus-nucleus collisions at SPS energies |0 have suggested that a reasonable 



description of the meson and baryon rapidity distributions can also be achieved on the basis 
of partonic degrees of freedom. However, in the latter calculations the extrapolation of the 
strong coupling constant to low has been overestimated as discovered recently by Bass 
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and Miiller This finding invalidates the detailed predictions and comparisons within the 
hybrid models VNI+UrQMD or VNI+HSD presented in Ref. H that have been tailored 
to describe the dynamics at RHIC or even LHC energies. 

We start with pp collisions at i/i = 200 GeV. The calculated results for the proton, 
TT"*" and rapidity distributions in the cms are shown in Fig. |T2] (upper part) for both 
models, which are denoted individually by the labels VNI and HSD in obvious notation. 
On the level of pp collisions we find only minor differences between the two kinetic models. 
The parton cascade shows a slightly higher amount of proton stopping as the HSD model 
(l.h.s.) and as a consequence a slightly higher production of 7t~^ and mesons (r.h.s.), 
because the energy taken from the relative motion of the leading baryons is converted to 
the production of mesons. It is presently unclear which of the two approaches will be closer 
to experiment; a proper description of pp data will be a necessary step before performing 
reliable extrapolations to nucleus-nucleus collisions. 

Inspite of this missing experimental information we directly step towards central collisions 
{b < 1.5 fm) for Au + Au at y/s = 200 GeV. The calculated results for the net proton (here 
p ~p), antiproton, 7r+ and rapidity distributions in the cms are shown in Fig. |T2| (lower 
part). In the HSD scenario essentially 'comover' scattering occurs with a low change of the 
meson rapidity distribution. Thus the meson rapidity distributions are roughly the same 
as for pp collisions. Also note that at midrapidity the net baryon density ^ Np — Np is 
practically zero, however, even at midrapidity at lot of baryons appear that are produced 
together with antibaryons. Thus also mesons (especially cc pairs) encounter a lot of baryons 
and antibaryons on their way to the continuum. Whereas the HSD approach predicts a 
vanishing net baryon density at midrapidity, other recent models - that combine high and 
low energy transport concepts - suggest a sizeable net proton density for Ucm ~ ||100|| . 

The amount of higher order hadronic rescattering processes at RHIC energies is depicted 
in Fig. |l^ (lower right part) as emerging from the HSD calculation, where the number of 
baryon-baryon [BB) and meson-baryon collisions {mB) is shown as a function of the invari- 
ant energy ^/s. We mention that quark-baryon and diquark-baryon collisions are counted 
here as niB or BB collisions, respectively. Apart from the initial small peak at ^/s = 200 
GeV a substantial amount of intermediate and low energy rescattering processes with max- 
ima at 2.5 GeV and 1.8 GeV are found, which essentially stand for flavor exchange processes, 
multiple pion production in mB and BB collisions as well as secondary strangeness produc- 
tion channels. For comparison the corresponding ^/s distributions are also displayed for 
bombarding energies of 2, 11 and 160 A-GeV, respectively, showing a dominance of low en- 
ergy BB and mB collisions. The latter reactions occur at energy scales where perturbative 
QCD is no longer applicable. This has to be kept in mind additionally when comparing to 
pp and pA reactions at y/s = 200 GeV. 

We return to the question of charmonium suppression at RHIC energies since it is ex- 
pected that one might probe increasing energy densities also with increasing centrality of 
the collision, where the latter can be correlated with the transverse energy Et produced in a 



collision event. As argued e.g. by Satz the survival factor S'j/* then should show steps 
as a function of Et due to the melting of first the Xc and then the J/\l/ in a QGP phase. As 
seen from Fig. 10 there are no pronounced steps in the Et dependence of J/\l/ suppression in 
the data for Pb -|- Pb at SPS energies according to the authors point of view; this situation 
might change at RHIC energies. 
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Using the 'late' comover model described in Section 5 we have calculated the J/^ survival 
factor Sj/is, as a function of the transverse energy Et in the cms rapidity window [-1,1] for 
Au + Au at a/s = 200 GeV on an event by event basis covering all impact parameters 
from 6 = to 13 fm. The resulting correlation of Sj/ii, with is shown in Fig. 14 and 
indicates a smooth decrease with centrality (or increasing Ex) reaching an average survival 
probability of ~ 0.1 for the most central events (cf. Fig. 11, r.h.s.). This result can be 
understood as follows: According to our calculations the net J/\l/ dissociation by mesons in 
central Au + Au collisions at the SPS is ~ 16% (cf. Fig. 11) while the rapidity distribution 
of negatively charged particles {h~) at midrapidity here is about 180. At the RHIC energy 
we get a corresponding h~ rapidity density at midrapidity of ~ 450 (cf. Fig. 12) which 
is higher by a factor of 2.5. Simply multiplying the J/\l/ meson absorption at the SPS of 
16% by the factor 2.5 we obtain about 40% for central collisions at RHIC energies, which 
together with ^ 52% of absorption on baryons gives a survival probability of 8%. The 
actual numerical results in Fig. 14 indicate that this simple estimate works quite well. On 
the other hand, if the h~ rapidity density is found to be lower (higher) experimentally, we 
expect corresponding changes in the J/\l/ suppression for central events if the 'late' comover 
absorption model holds true. This dependence might well be tested experimentally in the 
near future to possibly falsify the comover dissociation model. 

7 Summary 

In this work we have performed a systematic analysis of hadron production in central 
Au + Au collisions from SIS to RHIC energies within the HSD transport approach. We 
have concentrated here on the 'classical' signatures, i.e. strangeness and low mass dilepton 
enhancement as well as charmonium suppression. For all observables our calculations give a 
monotonic increase (for the ratio K^/ir^ and charmonium suppression) or decrease (for the 
low mass dilepton enhancement) with bombarding energy, respectively. So far, experimen- 
tal data are available only in a limited range of bombarding energies or at a single energy, 
respectively. We have pointed out that the relative maximum indicated by the experimental 
data in the /tt^ ratio at about 10 A-GeV (or higher?) is not reproduced within the 
transport approach that is based on quark, diquark, string and hadronic degrees of freedom. 
We speculate that at AGS energies this failure might be attributed to a restoration of chiral 
symmetry in a sufficiently large space-time volume (cf. Figs. ^ and ^). 

The enhancement of low mass dileptons in the range 0.3 GeV < Mg+e- < 0.6 GeV is most 
pronounced at lower bombarding energies of 2-5 A- GeV within our calculations since here 
the space-time volume for densities above 2 po is very large such that a majority of p-mesons 
probes the high density phase of the reaction (cf. Fig. ^). With increasing bombarding 
energy the average density - which a p-meson experiences - drops substantially such that 
high energy nucleus-nucleus collisions are not well suited for in-medium p spectroscopy. 

The suppression of charmonium (here J /'^) increases smoothly with bombarding energy 
and centrality of the reaction within the 'early' and 'late' comover absorption scenarios. 
Unfortunately, both scenarios cannot be distinguished by means of the excitation function 
since they give approximately the same survival probability Sj/ij, with bombarding energy 
(cf. Fig. 11). In the transport approach the smooth increase of charmonium absorption with 
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bombarding energy is easy to understand: a major fraction of J/\I/'s is anyhow dissociated 
by baryons which basically are of the same number at all energies considered here; only 
the relative collisional energy changes. The additional absorption by 'early' strings or 'late' 
hadrons increases smoothly with bombarding energy since the string and hadron density 
increases accordingly. At RHIC energies this additional suppression mechanism leads to a 
J/'^ suppression of about 90% in central Au + Au collisions even without employing an 
explicit formation of a QGP. We note, however, that the charmonium suppression shows a 
smooth dependence on the transverse energy Et (cf. Fig. 14); any gradual steps of Sj/^, 
with Et due to a melting of the Xc or the J/\l/ at higher energy density would indicate a 
new suppression mechanism which might be attributed to color screening in a QGP phase 
8l. 
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Figure 1: The nucleon scalar (Us) and negative vector potential {—Uo) as a function of the 
nuclear density p and relative momentum p of the nucleon with respect to the nuclear matter 
rest frame as implemented in the HSD transport approach [T^ . 
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Figure 2: The potential Usep (5) - as resulting from the nucleon scalar (Us) and vector 
potential (Uq) in Fig. 1 - as a function of the nuclear density p and relative momentum p of 
the nucleon with respect to the nuclear matter rest frame. 
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Figure 3: The scalar quark condensate < qq{x, 0,z;t) > for central Au + Au collisions at 6 
A -GeV divided by the vacuum condensate < qq>y such that the nonperturbative vacuum is 
characterized by a value of 1. The 2;-direction has been stretched by the Lorentz- factor ^cfn 
to compensate for Lorentz contraction, while negative numerical values for the condensate 
have been suppressed. 
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Figure 4: Same as Fig. 3 for central collisions of Au + Au at 20 A-GeV. 
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Figure 5: The meson (vr"*", rj, K^, K , (j) and J /'^) multiplicities from the HSD approach for 
central collisions of Au + Au from 200 A MeV to 21.5 A-TeV. 
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Figure 6: The j-n^ ratio at midrapidity from central Au + Au {Pb-\- Pb) collisions from 
1 A-GeV to 21.5 A-TeV. The open circles show the results from HSD for pp collisions while 
the open squares are obtained for central Au + Au reactions. The experimental data from 
Refs. [0, are displayed in terms of the full circles. 
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Figure 7: The differential p-meson distribution versus baryon density p/po (at the p creation 
point) for central collisions of Au + Au at 2, 5, 10, 20, 40, and 160 A- GeV from the 
HSD approach. The production channels involving two baryons or a meson and a baryon 
(denoted by BB p, nB p) are summed up in the solid histograms whereas the dashed 
histograms stand for the sum of the meson production channels (denoted by tttt p) which 
are dominated by the pion annihilation channel. 
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Figure 8: The differential dilepton multiplicity due+e- /dM for central collisions of Au + Au 
for 2, 3, 5, 10, 20, 50, 160, and 21500 A-GeV. Upper part: HSD calculations involving the 
'free' p-meson spectral function: lower part: HSD calculation involving the in-medium p 



spectral function from Rapp et al. |73 |. 
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Figure 9: The difFerential dilepton multiplicity due+e- / dM for central collisions of Au + Au 
at 2, 10, 50, and 160 A-GeV from the HSD calculations involving the 'free' p-meson spectral 
function (solid lines) and the in-medium p spectral function from Rapp et al. |^ ^ (dashed 
lines) for comparison. 
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Figure 10: The J/\l/ suppression (in terms of the fi^fi^ decay branch relative to the Drell- 
Yan background from 2.9 - 4.5 GeV invariant mass) as a function of the transverse energy 
Et in Pb + Pb colhsions at 160 A-GeV. The sohd hne stands for the HSD result within 
the 'late' comover absorption scenario from Ref. [|85| while the dashed line results from the 



'early' string absorption scenario from Ref. involving a transverse string radius 
0.2 fm. The full dots stand for the NA50 data from 1995 



33 1, the full squares for the 1996 



data |7^, the open triangles for the 1996 data with minimum bias while the open circles 
represent the 1998 data |^4[. Note that the 1995 data have been rescaled in Et as compared 
the same rescaling has been adopted to the calculations from Refs. [p5l 



to Ref. []63 



which had been compared to the data from . 
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Figure 11: The J/"^ survival factor 5*"^/* (in terms of the jJi'^jJ.^ decay branch relative to the 
Drell-Yan background from 2.9 - 4.5 GeV invariant mass normalized to the same quantity 
for pd reactions) as a function of the bombarding energy in central Au + Au collisions from 
20 to 160 A-GeV (l.h.s.) and from 160 A-GeV to 21.5 A-TeV (r.h.s.). The solid line (full 
dots) stands for the HSD result within the 'late' comover absorption scenario from Ref. 
while the solid line (open circles) results from the 'early' string absorption scenario from 
Ref. involving a transverse string radius = 0.22 fm in order to match both absorption 
scenarios at 160 A-GeV. The dashed lines stand for the relative fraction of J/\E' dissociations 
with baryons while the dotted lines stand for the 'early' comover string dissociation (open 
circles) and 'late' comover meson dissociation (full dots), respectively. The total suppression 
factors (full lines) are practically identical for both scenarios from 20 to 160 A- GeV (l.h.s.). 
This also holds for the energy range from 160 A-GeV to 21.5 A-TeV within the numerical 
accuracy achieved. 
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Figure 12: The rapidity distribution of protons (upper left part) from pp collisions at = 
200 GeV from the HSD approach (solid histogram) in comparison to the prediction from the 
parton cascade VNI ||9^ (dashed histogram); the upper right part shows the same comparison 
for the 7r+ and rapidity distributions, respectively. The lower part of the figure displays 
the HSD predictions for central {b < 1.5 fm) Au + Au at y/s = 200 GeV per nucleon; (l.h.s.) 
net proton {p — p) and p rapidity distribution, (r.h.s.) 7i~^ and rapidity distributions. 
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Figure 13: The number of baryon-baryon (solid histograms) and meson-baryon colhsions 
(dashed histograms) as a function of the invariant energy ^/s for central {b < 1.5 fm) Au + Au 
collisions at bombarding energies of 2, 11, 160 A-GeV and 21 A-TeV, respectively, from the 
HSD model. The arrows indicate the string thresholds for mB and BB collisions of 2.3 GeV 
and 2.6 GeV, respectively. 
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Figure 14: The J/\l/ survival factor 5*7/^, as a function of the transverse energy Et in the 
rapidity interval (-1 < ycm < 1) in Au + Au collisions at ^/s = 200 GeV in the 'late' comover 
model. The solid line represents the result for J/\l/ dissociation on nucleons and mesons, 
whereas the dashed line and the dotted line correspond to the absorption on baryons and 
mesons, respectively. The error bars in the figure are due to statistics only. 
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